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Abstract. We investigate some modifications to the 
Press & Schechter (1974) (PS) prescription resulting from 
shear and tidal effects. These modifications rely on more 
realistic treatments of the collapse process than the stan- 
dard approach based on the spherical model. 

First, we show that the mass function result- 
ing from a new approximate Lagrangian dynamic 
(Audit & Alimi 1996), contains more objects at high 
mass, than the classical PS mass function and is well fitted 
by a PS-like function with a threshold density of 6c — 1.4. 
However, such a Lagrangian description can underesti- 
mate the epoch of structure formation since it defines it 
as the collapse of the first principal axis. We therefore 
suggest some analytical prescriptions, for computing the 
collapse time along the second and third principal axes, 
and we deduce the corresponding mass functions. The col- 
lapse along the third axis is delayed by the shear and the 
number of objects of high mass then decreases. 

Finally, we show that the shear also strongly affects 
the formation of low-mass halos. This dynamical effect 
implies a modification of the low-mass slope of the mass 
function and allows the reproduction of the observed lu- 
minosity function of field galaxies. In a companion paper, 
we present results of numerical simulations which com- 
plete this work. 
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1. INTRODUCTION. 

One of the simplest statistical indicators which character- 
izes the mass distribution in the universe is the mass func- 
tion of cosmic structures. Press and Schechter (1974, here- 
after PS), have proposed a formalism to evaluate this func- 
tion at any time. For a given hierarchical scenario of struc- 
ture formation, as for instance CDM, the mass function 
is obtained using only the power spectrum of the initial 



density fluctuations. This formalism has been extensively 
tested with the help of numerical simulations. The first 
comparison between the mass function of dark-matter ha- 
los obtained in numerical simulations and that predicted 
in the framework of the PS approach have shown a satis- 
factory agreement (Efstatiou et al. 1988). This agreement 
was later improved by imposing some ad-hoc modifica- 
tions to the initial PS formalism (Efstatiou & Rees 1988, 
Calberg & Couchman 1989, Lacey & Cole 1994). For ex- 
ample, it was shown in these studies that the concordance 
between the numerical and the analytical predictions was 
improved by lowering the initial density threshold at which 
a spherical homogeneous region collapses. However, the 
theoretical reasons for which a quantity resulting from 
very non-linear dynamics, such as the mass function, can 
be computed from this formalism have never been clar- 
ified. This question takes on a greater importance be- 
cause the PS formalism which is a powerful tool, is widely 
used in phenomenological analytical cosmology. It is, for 
example, very often used to put constraints on the cos- 
mological parameters using either the luminosity function 
of galaxies (White & Frenk 1991, Blanchard et al. 1992), 
the temperature 
distribution of X-ray clusters (Oukbir & Blanchard 1992, 
Viana & Liddle 1996) or the Sunayev-Zel'dovich effect 
(Cole & Kaiser 1988, Barbosa et al. 1996, Eke et al. 1996 
and references therein). 

The mass function of cosmic structures obtained in 
the framework of the PS formalism has a statistical and 
dynamical origin. Starting with the statistical properties 
of the initial density field, PS identify the fraction of mass 
embedded in collapsed objects of mass greater than a given 
mass M, with the fraction of volume filled with regions 
of space satisfying a given dynamical criterion. Usually 
the criterion is simply issued from the spherical collapse 
which selects all the points whose initial density is above a 
certain threshold. It is highly probable that such a simple 
model is unable to describe the complex features occurring 
during a fully non-linear dynamical process. Therefore, 
one can imagine that only a fraction of regions with a 
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given density contrast exceeding the threshold eventually 
collapses. This "fuzzy" density criterion can be viewed as a 
statistical tool (Blanchard et al. 1992), but here, in order 
to have a better understanding of structure formation, we 
intend to derive it from the gravitational dynamics. 

Therefore, in this paper, we study the influence of non- 
linear dynamics on the mass function, highlighting the role 
played by shear and tidal effects on the actual shape of the 
mass function. We present analytical results which exhibit 
generic features of PS-like mass functions and show the 
relationship between the dynamics and the mass function 
through what we call a selection function. Our study con- 
cerns both non-spherical Lagrangian dynamical models, 
and some dynamical "ansatz" that we introduce to de- 
scribe the structure formation beyond the usual validity 
range of Lagrangian approximations. In a companion pa- 
per we compare our analytical predictions with numerical 
results. 

2. Mass Functions from Lagrangian Dynamics 

2.1. PS formalism from a Lagrangian point of view. 

The analytical method which has been most successful up 
to now to compute the mass function is the one intro- 
duced by PS. In this paper we use such a formalism but 
we intend to give our own interpretation of this approach 
from a Lagrangian point of view (see also Monaco 1995). 
This gives us a general framework for computing the mass 
function. 

The idea of PS is to compute the number density of 

cosmological objects with mass M at a given time, by 
identifying the fraction of mass contained in these objects 
with the mass embedded in the space filled by points with 
a density which reaches infinity at the same time. When 
considering an object of mass M, the initial density field 
is smoothed on a scale R with M oc R^.To give a concrete 
expression for this idea one needs a dynamical prescription 
that allows to follow the evolution of each point of the 
initial density field and to know its statistical properties. 

Throughout this paper we work in a Friedman- 
Robertson- Walker universe with a critical density param- 
eter, Sl = 1, and initial Gaussian random fluctuations. At 
some early epoch, corresponding to an expansion factor 
ai and a redshift Zi, the density fluctuations are entirely 
described by their power spectrum. We flrst multiply this 
spectrum with a low-pass filter of characteristic scale R 
that cuts the power above a frequency ~ 2tt/R. The choice 
of fllter is, a priori, arbitrary (Lacey & Cole 1994). We re- 
strict ourselves to the Top-Hat window function in real 
space, because its application has an intuitive physical in- 
terpretation. It allows to consider homogeneous and spher- 
ical regions of space. We then assume that the resulting 
density fleld is correctly described by independent fluid 
elements. Each fluid element is now initially described 
by some dynamical parameters, which are obtained by 



smoothing the corresponding fields to the scale R. Given 

these initial conditions, we evolve each fluid element sep- 
arately until it collapses (i.e. until the density of the fluid 
element becomes infinite). The evolution has to be fol- 
lowed in a Lagrangian way in order to relate the collapse 
epoch of the fluid element to its initial conditions. In the 
standard PS approach, these fluid elements evolve accord- 
ing to the homogeneous spherical model. In this case, the 
collapse epoch is entirely determined by the initial density 
contrast of the fluid element. However, it has been shown 
by several authors (e.g. Van de Weygaert & Babul 1994, 
Audit & Alimi 1996) that, during the gravitational col- 
lapse, the shear and the tidal forces are very important. 
Consequently, more complex, non-spherical, gravitational 
dynamics have to be considered in order to take into 
consideration these physical effects which are totally dis- 
carded in the spherical model. 

For Lagrangian non-spherical dynamics, the collapse 
time Oc of a fluid element can depend on several initial 
smoothed parameters x^, j = l,n. The fraction of mass 
embedded at present (a = 1) in collapsed structures with 
a mass greater or equal to M is then given by 

F(> M) = j Pr (x^) s (x^) dx^--- rfa;" (1) 

where s = 1 if Cc (x^) < 1 and otherwise. The selection 
function s makes a direct link between the Lagrangian dy- 
namics of the fluid elements and the mass function. The 
integral of s, with respect to all Xj parameters except the 
initial density contrast, can be interpreted as a statistical 
quantity. This quantity gives the probability that a given 
region of space becomes part of a bound structure accord- 
ing to its initial density contrast (Blanchard et al. 1992). 
The right-hand side of equation (1) depends on M through 
the smoothing scale R. Pr is the probability density func- 
tion of the initial smoothed parameters, it represents the 
fraction of volume fllled by the regions of space having 
initial parameters {x^ , j = 1 , n) . 

In general, the above definition runs into the so called 
"renormalization problem", i.e. F{> 0) = Fq 7^ 1. 
This problem can be alleviated if one carefully consid- 
ers the fate of regions which have not been selected on 
a scale i?,, but can be selected with a non-zero proba- 
bility on another scale R' > R, and therefore belong to 
objects with a mass greater than M. Excursion set the- 
ory (Peacock & Heavens 1990, Bond et al. 1991a, 1991b) 
solves this renormalization problem for the particular case 
of the sharp-k filter, but the physical interpretation of this 
filter remains unclear. In this work we do not investigate 
this renormalization problem. We always assume that this 
renormalization is justified in any case. We therefore di- 
vide the mass fraction obtained in equation (1) by the 
constant Fq. The mass function of collapsed objects, which 
gives the number density of object with mass M, is then 
defined as 
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$(M) = - 



po d (F{> M) 
MdM 



(2) 



In order to determine the mass function we therefore 
have to cvahiatc the statistical properties of the initial 
parameters and the collapse epoch associated with these 
parameters. In particular, we evaluate the collapse epoch 
using three different Lagrangian dynamics: the spherical 
model; the Zel'dovich approximation and a more realistic 
local Lagrangian model that we introduced in a former pa- 
per (Audit & Alimi 1996). These three dynamical models 
state that the fluid element belongs to a collapsed ob- 
ject when its density becomes infinite. For the two latter 
cases, this first singularity corresponds to the collapse of 
the fluid element along its first principal axis and then 
favors the formation of sheet-like objects. One can then 
ask whether such objects agree with the idea of having 
a dense, virialized "blob" in the density field. Unfortu- 
nately, Lagrangian theory is unable to bypass this first 
singularity. We therefore investigate two other models in 
section 3 which roughly describe the collapse of the fluid 
element along its second and third axis. This comparative 
analysis of the different mass functions we obtain for each 
case allows one to study the effect of the dynamics, and in 
particular, the influence of the shear and the tide on the 
formation of structures in a PS formalism. 

2.2. Computation of the Collapse Epoch 

In the simplest cases, only one initial parameter is needed 
to determine the collapse epoch of a fluid element. This 
is the case for the spherical model, where the initial den- 
sity contrast entirely determines the evolution of a fluid 
element, and for the Zel'dovich approximation where the 
smallest eigenvalue of the deformation tensor gives the 
collapse epoch. However, these dynamics are rather sim- 
ple. It does not seem very realistic to account for the va- 
riety of physical processes occurring during the collapse 
with a single initial parameter. Therefore, we also use 
a local Lagrangian dynamic which takes into account, 
to some extend, the effect of the shear and the tide 
(Audit & Ahmi 1996), and then determines the collapse 
epoch in terms of the three eigenvalues of the deformation 
tensor. 

The collapse epoch (written here as the collapse red- 
shift Zc) of a fluid element evolving according to the spher- 
ical model is given by the simple formula 



S_ 

5c 



(3) 



where 5c = t,Vt) - 1-686. We use for aU fields its 
initial value linearly extrapolated to = (for example: 
5 = {1+Zi)5{zi)). Within the spherical model , the collapse 



epoch depends only on the initial density contrast and 
only fluid elements which are initially over-dense, collapse. 

In the framework of the Zel'dovich approximation 
(Zeldovich 1970) the density contrast of a fluid element 
evolves according to 



1 + 5{t) = 



1 



(1 + a(i)Ai)(l + a{t)X2){l + a{t)X3) 



(4) 



where Ai, A2 and A3 are the three initial eigenvalues of 
the deformation tensor arranged in decreasing order. The 
collapse epoch of the flrst principal axis is reached when 
the denominator in equation (4) vanishes 



1 + = -A3 = g - 0-3 



(5) 



Once again S is the initial density contrast of the fluid 
element and are its initial shear eigenvalues. In the case 
of zero shear, the collapse epoch is given by: 1 + Zc = 5/3. 
This means that the exact solution of the spherical case 
is not recovered, as can be expected for this first order 
dynamical prescription. For a fixed density contrast, the 
shear accelerates the collapse (0-3 < 0) and equation (5) 
can be considered as a first simple approximation for non- 
spherical dynamics. 

The Zel'dovich approximation illustrates in a simple 
way that the collapse epoch cannot be considered as a 
function of the density contrast alone. However, comput- 
ing the collapse epoch in terms of A3 alone is not satis- 
factory either. Therefore, we present in this paragraph a 
derivation of the collapse time which uses a Lagrangian 
dynamic based on the deformation tensor. Audit & Alimi 
(1996), hereafter AA, have proposed a local Lagrangian 
dynamic to determine the evolution of the eigenvalues of 
the deformation tensor. Knowing these eigenvalues and 
their time derivatives is enough to compute all the other 
dynamical quantities. The dynamical prescription of AA 
is well suited to the computation of the collapse epoch. It 
accounts well for the effect of the shear on the collapsing 
fluid element. It is exact for spherical, cylindrical and pla- 
nar geometry and very accurate over a wide range of ini- 
tial conditions. Similar to the Zel'dovich approximation, 
it predicts that under-dense regions can collapse due to 
the effect of the shear. 

In terms of the eigenvalues of the deformation tensor, 
the density contrast is defined as 



■5{r) 



1 



(1 + Ai(r))(l + A2(T))(1 + A3(r)) 



The evolution equation for each of the Xi{T) is given 
in the framework of the AA approximation by 



l + A,/2 + Afc/2 + AjAfc/2 \ 6A£ 

1 + Xj + Xk + XjXk ) ^ ^ 
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where r is the conformal time defined by: dr = dt/a^ and 
(i,j,k) is a circular permutation of (1,2,3). These three 
equations are highly non linear and have been integrated 
numerically. 

The initial eigenvalues of the shear and deformation 
tensors are related by: Xi = ai — S/3. The two triplets 
(Ai, A2, A3) and (i5, a2, u^) are therefore totally equivalent 
(tTi is obtained using the relation: cti + (T2 + (73 = 0). 
In the following, we prefer to use the triplet {6,a2,cr^), 
which allows an immediate physical interpretation. The 
eigenvalues of the shear and deformation tensor are always 
arranged in decreasing order: cti > (T2 > (T3 and Ai > A2 > 
A3. 

Moreover, in order to integrate equation (6) we only 
keep the linear growing mode in the initial conditions. 
This restriction and the invariance of equation (6) under 
the transformation: t — > ar ensures the following scaling 
law for the collapse epoch 



two full lines gives an idea of the influence of (T2 on the dy- 
namics and of the importance of the coupling between the 
different directions of collapse of the fluid element. Note 
however that the error one could make by neglecting the 
influence of (T2 to compute the collapse epoch of the first 
principal axis is of the order of 20%. 



Oc ((5, 0-2 , 0-3 ) = offlc {a6, aa-2. , acr3 ) 



(7) 



Such a property of the collapse time is also present for 
the spherical model and the Zel'dovich approximation. It- 
allows to reduce by one dimension the parameter space 
that we have to explore to determine the collapse epoch 
for any initial condition. It is therefore sufficient to know 
the function T+(o-2,o'3) (resp. (0-2, 173)) representing 
the inverse of the collapse epoch (a^^ = 1 + Zc) in the 
(f2,C3) plane for i5 = -|-1 (resp. 6 = —1). Then, we can 
deduce the collapse epoch for any set of initial conditions 
using the scahng law of equation (7). In figure (1), we 
plot the iso-contours of the functions T+ and T~. From 
the previously given relations between the ai it is easy to 
deduce that (J3 < a2 < — (T3/2. This explains why only a 
portion of the (0-2, 0-3) plane is covered. For each of these 
plots, a thousand points were computed by numerically 
solving equation (6). We present an analytical fit (accurate 
to within 1%) of these two surfaces in Appendix A. For 
(5 > 0, the fluid element always collapses in a finite time. 
For i5 < 0, a finite area of shear values is singular (i.e. 
the collapse time is infinite). This area corresponds to a 
fiuid element whose initial deformation tensor has three 
positive eigenvalues (all fiuid elements with at least one 
negative eigenvalue collapse). 

The two lower panels of figure [1] show the inverse 
of the collapse time as a function of (T3 for S — +1 and 
5 — —1. The long-dashed fine corresponds to the spheri- 
cal model for which the collapse time is independent of . 
The small-dashed line corresponds to the Zel'dovich ap- 
proximation. The full lines represent the inverse collapse 
time of the AA dynamics. The upper line corresponds to 
the case where (T2 = 0-3 (the fastest collapse at given 0-3 
and S) and the lower line to (T2 = —^3/2 (the slowest col- 
lapse at given 0-3 and S). For 6 < the collapse time be- 
comes infinite as (73 reaches 6/3. The spread between the 




Fig. 1. Upper panels: Iso-contours of the redshift of collapse 
in the (0-3, (J2) plane for S = +1 and 5 = —1. The shaded 
area corresponds to fluid elements which never collapse. Lower 
panels: The collapse redshift (or the inverse of the collapse 
epoch) as a function of the smallest eigenvalue of the shear. 
The solid curves correspond to the two boundary lines of the 
triangles plotted in the upper graphs. The small-dashed lines 
correspond to the Zel'dovich approximation. For S = +1 the 
spherical model is shown as a long-dashed line. 



2.3. Statistics of the Initial Parameters 

In the case of initial Gaussian density fiuctuations and 
assuming an initial linear growth mode only, Doroskevitch 
(1970) derived the probability function for the eigenvalues 
of the deformation tensor, namely 



^'(Ai, A2, A3 



675V5 

e A 

87rA6 

x(Ai - A2)(Ai 



1532 
2A^ 



A3)(A2-A3) (8) 
where A is the r.m.s density contrast smoothed at scale R 



A^{R) = J 4Trk^P{k)W^{kR)dk 



(9) 
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and si = Ai + A2 + A3 and S2 = A1A2 + A2A3 + A1A3. The 

A's arc correlated variables, S and the doublet ((72, (73) are 
independent. The probability distribution for these new 
variables can then be deduced easily from equation (8) 



/CT2 CT3\ 

VA' a) 



(10) 



where the marginal probability functions for the rescaled 
variables v = bj A and \Xi = Oij A are given by 



1 



-1/^2 



27r 



(11) 



and 



^ , / Ap-7-5(/i^+M2M+/:*l) 

4 V27r 

X (/i2 + 2/U3)(/i3 + 2/X2)(M2 - Ms) 



We compute also the marginal probability distribution 
functions of /ii e [0, -|-oo[ and /U3 e] — 00, 0] 



(12) 



and of H2 s] — 00, +00 [ 



-45,*^ /2 



(13) 



Equation (1) can be now re- written in terms of u, ^2, 
Us and A. 



F{> M) = jp^Pi_,^^i^^s{iyA, 112A, ij,3A)diydti2dn3 



(14) 



The only quantity which then depends on M is A. 
Consequently, the mass function can be expressed as 



$(M) = -^^J- 

/ds 
P'^Pn2,n3-g^{'^^^ M2 A, iJ,3A)dudfi2dfi3 (15) 

where the integral is a function of A and depends only on 
the chosen dynamics which gives its specific form to s. It is 
therefore possible to define a universal mass fimction for 
each dynamic which is independent of the initial power 
spectrum of the density fluctuation. The mass function 
can be written, for any spectrum, in the form 



$(M) 



Po dA 
' MdM 



$(A) 



(16) 



The spectrum dependence of the mass function is to- 
tally contained in the dA/dM factor. We call ^'(A), which 
depends only on the chosen dynamics (i.e. on the func- 
tion 1 -I- Zc), the universal mass function. This demon- 
strates that the mass function resulting from a PS for- 
malism admits, in addition to the well know time self- 
similarity, a "scaling" behavior in spectrum (see also 
Lacey & Cole 1994). This property holds only in a critical 
universe, independently of the dynamics as long as only 
the growing mode is kept during the linear regime. These 
two "scaling properties" of the mass function are general 
characteristics of a PS-like mass function. Therefore, they 
can be used as test for the validity of a PS approach to 
the mass function. 

2.4- Mass Functions 

Equation (1) can now be re- written as follows 



1 p + <x> 

F{> M) = -= / e-'' l'^S{yA,A)dv 
V27r J_(x) 



(17) 



For the spherical model, S is identical to the selection 
function s which depends then only on the initial density 
contrast 



s(i/A) 



if !/A < (Sc 

1 otherwise 



(18) 



For the Zel'dovich and AA approximations, the se- 
lection function is not a function of the density contrast 
alone. In equation (17), S is then the selection function, 
integrated over all parameters (/X3 for Zel'dovich approx- 
imation and \X2 and for AA approximation), except v. 
For the AA dynamical model, S is explicitly given by 



S{vA,A) = f_^\Xl^ 5(i.A,M2A,M3A) 



xP„ 



jdii2dfi3 



(19) 



It represents the probability that a fluid element with a 
given initial density contrast, collapses in the age of the 
Universe. In figure (2), we plot this quantity and its deriva- 
tive with respect to A for A = 0.5, 1 and 2. A small 
A value corresponds to high mass objects, and a large 
A value to small mass objects. One clearly sees that for 
A = 2, negative values of the density contrast are selected, 
showing that shear effects can actually lead to the collapse 
of under-dense regions. For A = 0.5, the selection func- 
tions become steeper and converge towards an Heavyside 
function. For the Zel'dovich model, the selection function 
in the limit A — > leads to a Heavyside step function 
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at 5 = 3, illustrating that this model is only a first or- 
der approximation. For the AA model, the selection func- 
tion converges in the limit A ^ towards the Heavyside 
function at S = 6c- The derivative of the selection func- 
tion dS/dA peaks at a value lower than 5c, and this peak 
value tends towards (5c as A ^ 0. Therefore, in the AA ap- 
proximation, where objects are defined by the collapse of 
their first principal axis, the mean initial density contrast 
is lower than the standard spherical threshold 6c- 
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Fig. 2. Selection functions and their derivatives (see text for 
definitions) obtained at A=2 (dotted lino), A=l (solid line) 
and A=0.5 (dashed line) for the Zel'dovich model (two lower 
graphs) and for the A A model (two upper graphs). 

We deduce now the universal mass function 

For the spherical dynamics Fq = 0.5 since only 50% of 

the regions of space collapse (those with a positive den- 
sity contrast). For the Zel'dovich and AA approximations, 
Fo ~ 0.92 as all regions with a negative A3 collapse. In fig- 
ure (3) we plot $(A) for the three different dynamics. The 
bottom figure is shown with a linear scale to allow an easy 
examination of the maximum amplitude of the mass func- 
tions. The top one, shown on a logarithmic scale, is better 
adapted to observing the small mass (high A) slope and 
the high mass (small A) cut-off. The mass function com- 
puted with our dynamical model (full line) is very well 



fitted by a standard PS mass function with a critical den- 
sity parameter (5c = 1.4 (short-dashed line). This is easily 
explained by the behavior of the selection function previ- 
ously discussed, whose effective threshold is lower than the 
spherical threshold 6c- On the contrary, for the Zel'dovich 
approximation, the resulting mass function cannot be ap- 
proximated by a PS function with another value of ^c- The 
high-mass cut-off is however located at lower mass, as can 
be expected from the selection function which chooses re- 
gions with an effective threshold higher than the standard 
value 5c- All mass functions have an exponential cut-off 
at high mass and a power law behavior at small mass. 
However, this doesn't mean that the resulting mass func- 
tion can always be fitted by the PS formula, with 6c as fit 
parameter. 




Fig. 3. Universal multiplicity functions on linear and logarith- 
mic scales. The solid line corresponds to our dynamical model 
(Audit & Alimi 1996), the dashed line to the PS mass function 
with 5c = 1.4, the dotted line to the standard PS mass func- 
tion (5c — 1.686) and the long-dashed line to the Zel'dovich 
approximation. 

Moreover, we see that the standard {5c — 1.686) PS 
mass function (dotted line) and our mass function have 
roughly the same slope for small mass objects, although 
we predict more high mass objects than the standard PS. 
For example, if the power spectrum is normalised with the 
condition A(i? = 8Mpc.h~^) = 1 and assuming assume 
a critical Univers with h = 0.5, we find that our mass 
function contains 30% more structures of mass « IO^^Mq, 
which is the typical mass of galaxy clusters. 
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3. Beyond LagreingiEin Models 

In the previous section we have computed the mass func- 
tion in the framework of the PS formahsm using different 
Lagrangian dynamical models to determine the collapse 
epoch of a fluid element. In all cases, the collapse was 
defined as the epoch when the density contrast becomes 
infinite. It is then impossible to follow, with Lagrangian 
dynamics, the evolution of a fluid element beyond this flrst 
singularity. Moreover, if one considers that the fluid ele- 
ments describe extended regions of space, this first singu- 
larity corresponds to a stage where the region has totally 
collapsed along a single direction and reached a sheet-like 
geometry. The resulting "pancakes" may, however, not be 
the dense virialized halos that one would expect as well 
defined objects in the density field. It is therefore inter- 
esting to investigate what happens when one waits for 
the collapse of the fluid element along its second or third 
axis which corresponds to filamentary or quasi-spherical 
objects respectively. The mass function deduced from the 
collapse of the third axis will be presented in greater de- 
tail, as it is, in our opinion, best suited for dense virialized 
objects in the density field. 

In section 2.3 and figure (1), we have seen that the dis- 
persion of the collapse epoch of the first axis, in the AA 
approximation, due to the shear on the second axis was of 
the order of 20%. This weak dependence suggests that, as 
in the Zel'dovich approximation, the different directions 
of the fluid element can be considered as collapsing in- 
dependently. This approximation allows to artiflcially go 
beyond the first singularity and to propose an ansatz for 
the collapse epoch of the second and third axis in order. 

We require that the collapse epoch of a fiuid elements 
along any of its principal axis satisfies the following prop- 
erties: 

— The collapse epoch along a given axis is a function of 
the corresponding initial shear eigenvalue and of the 
initial density contrast. 

— The scaling properties (eq. [7]) should be satisfied. 

— The spherical collapse epoch should be recovered if the 
initial eigenvalue of the shear along this direction is 
zero. 

— The collapse epoch should be a monotonically increas- 
ing function of the shear. A positive (resp. negative) 
shear eigenvalue should slow down (resp. accelerate) 
the collapse relative to the corresponding spherical 
model. 

These properties mean that the following formulae hold 

ac{S,(Ti) = ^ac{l,ai/S) = ^dc{ai/d) (21) 

where he is a monotonous increasing function and adO) = 
5c- Since oLc is monotonous, it is also bijective and has an 
inverse, ctc, which gives the shear corresponding to a given 
collapse epoch in the 5 = +1 plane (note that there are, 



as in the previous section, two functions cLc, one for the 
positive and one for the negative density contrasts). 

3.1. Mass Functions of Sheets, Filaments or Knots 

In order to push our investigations further we now pro- 
pose a simple ansatz for the collapse time along each axis. 
We do not have an underlying dynamical model to justify 
for the following approximation of the collapse epoch of 
the different axes. These models are a phcnomcnological 
approach to studying qualitatively the infiuence of consid- 
ering the collapse along a given axis. 

For (5 > and (7^ < as well as for (5 < we approxi- 
mate by linear functions the results given in the previous 
section 

l + Zc = S/Sc -aui if (5 > and (Tj < (22) 
\^Zc = a{5/Z -Ui) if J < (23) 

With this model we recover, for i5 > 0, the spherical col- 
lapse epoch for cTj = and for (5 < 0. The collapse time 
becomes infinite for ai = 5/3. The slope a is a free pa- 
rameter, it has to satisfy 0.8 < a < 1 in order to recover 
approximately the right collapse time for the first axis (fig. 
4). 




-1 pi -1.5 -1 -0.5 



Fig. 4. The collapse redshift as a function of the eigenvalue 
of the shear. The ansatz defined by Eqs. (22) and (23) are 
shown by the full-line for the part deduced from the linear fit 
of the results of previous section. The dashed-line represents 
the extrapolation in the 5 > Q,(Ji > Q region (a = 0.8, e = 1) 
(Eq. 24). As a reminder, the results of the previous section are 
shown as a dotted-line. 

To compute the mass function for any axis, we also 
need to know the collapse epoch for 5 > and ai > 0. 
Unfortunately, in this case we do not have any dynamical 
model to determine the function dc- However, to study 
the qualitative behavior of the mass function we choose a 
simple ansatz for the collapse time given by 

l + Zc = Ul-^^r (24) 

Oc CTcO 
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The spherical collapse time is once again found for vanish- 
ing shear The shear totally prevents the fluid element from 
collapsing when ai > ac6. Note that for large values of ac 
we find the spherical model everywhere, e is a positive, 
free parameter which determines the global shape of the 
function he- This parameterization is somewhat arbitrary 
but satisfies the general requirements presented above. We 
now compute the selection function, its derivative and the 
universal mass function for objects which collapse along 
the first, second or third axis. 

3.1.1. First axis collapse 
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Fig. 5. Derivative of the selection function for A = 0.5, 1 and 2 
(long-dashed, solid and dotted lines respectively) and for first, 
second and third axis collapse (fig. a, b and c respectively). 
Figure d shows the selection function for A = 1 for the col- 
lapse of the first, second and third axes (long-dashed, solid and 
dotted lines respectively). In all figures the vertical solid line 
indicates the position 5 = 5c. 



This case has been studied in detail in the previous 
sections. We want to show here that the results are not 
strongly affected in neglecting the influence of a2 ■ Since 0-3 
is always negative we use formulae (22) and (23) for the 
collapse epoch. We have plotted in figure (5) the selection 
function and its derivative for the case a = 0.8 and for 
different values of A. The difference with the results of 
the previous sections is always less than a few percent. 
As was already mentioned before, the selection function 
is equal to unity for 5 > 5c and continuously decreases 
towards zero when ^ — > —00. Its derivative sharpens and 



peaks at a value which goes towards 5c (while remaining 
smaller) as A decreases. The resulting mass function is 
plotted on figure (6) and always differs by less than 10% 
to the one obtained in the previous section. 

3.1.2. Second Axis Collapse 

The collapse epoch of the second axis is a function of 5 and 
(72. In the last section, (T3 varies between —00 and 0. Now, 
(72 belongs to the domain ] — oo,+oo[. Therefore we use 
for the collapse epoch, the formulae (22), (23) and (24). 
We have plotted in figure (5) the selection function and 
its derivative for the case a = 0.8, e — 1 and ctc — l/aSc 
(we have just extrapolated formula (22) for the domain 
a2 > 0). In this case the selection function increases from 
zero to one when S goes from —00 to -l-oo and is equal to 
one-half for S — Sc (all the points with /Z2 < (50%) which 
collapse faster than the spherical model are selected). As 
before, a fraction of the fluid elements with 5 < 5c are 
selected, but contrary to what happened in the previous 
case, a fraction of the fluid elements with 6 > 5c are unable 
to collapse because of shear effects. The derivative of the 
selection function peaks at around 5c (exactly at 5c for our 
model) and has a width which tends towards zero in the 
limits of vanishing A. 

The resulting mass function is plotted in figure (6). 
It is now quite close to the standard PS mass function, 
which means that the number of objects with A < 1 has 
considerably decreased compared to the first axis collapse. 
This fact is not very surprising because in this case, the 
shear can either slow down the collapse for a2 > or 
accelerate it for 0-2 < 0. These two effects roughly cancel 
each other, and one obtains a quasi-spherical model. 

3.1.3. Third Axis Collapse 

We now present in greater detail the mass function ob- 
tained for objects that collapse along their third axis. 
In this case the collapse epoch is a function of 5 and 
(Ti e [0,+oo[. It is possible to derive analytically the se- 
lection function, its derivative and the corresponding mass 
function for a given dc- After some algebra, we obtain 



55„(i/A,A) = / P^Jn, 
Jo 



(25) 



Taking the derivative of this equation with respect to A 
yields 



{uA, A) = P^, {lyaduA)) v'^dc{vA) (26) 



where dc symbolizes the derivative of the function dc ■ Fi- 
nally, using the variable x defined by a; = Uciyt^^ the 
universal mass function can be expressed as 
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and for the variance 
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Fig. 6. Mass function obtained for first, second and third axes 
coUapse (dotted, dashed and long-dashed lines). The solid line 
is the classical {5c = 1.686) Press and Schechter mass function. 



where ac s]0, -|-oo] is the value of x at which Uc becomes 
infinite. A similar calculation can be carried out for all 
axes. The above formula is a general relation between the 
mass function and the collapse time defined by the under- 
lying dynamical model. 

For the third axis, we assume that the collapse epoch 
is now given by equation (24) and we use the same param- 
eters as for the second axis. The selection function (Fig. 
[5]) tends to zero as S ^ 6c and to unity as (5 ^ -|-oo. 
This reflects the fact that a fraction of the fluid elements 
with 6 > Sc are unable to collapse as the shear inhibits 
the collapse. 

The derivative of the selection function behaves sym- 
metrically about the line S = Sc compared to the one ob- 
tained for the first axis. Its derivative sharpens and peaks 
at a value which tends towards 6c (while remaining larger) 
as A decreases. In the class of models with e = 1, it is pos- 
sible to derive an analytical expression for these qualita- 
tive features. We compute the mean (which is close to the 
peak value) and the variance (which corresponds to the 
width) of the function dS/dA. We obtain for the mean 

3 r5~A A 

S = 6c + ^\hi ^<5c + 0.54— (28) 

5 V 27r (Tc (Tc 



^ 1 /297r-81 A ^^„A 

S=-a/ — — -0.19— (29) 

SV IOtT (Tc (Tc 

In these last two equations the width and the peak posi- 
tion are linear functions of A. Note also that as A — > 0, 
we get 6 = Sc and E = 0, as in the spherical model. 

The mass function can also be derived analytically. 
Using formula (27), we obtain 



$(A) = 




The mass function now has two "dynamical parame- 
ters", namely 6c and cTc. The first one describes the influ- 
ence of the initial density contrast on the collapse epoch 
and the second one gives the degree of inhibition of the 
collapse by the shear and the tidal field. The spherical 
mass function is recovered in the limits CTc — > +oo and for 
small values of (Tc most of the fluid elements are unable to 
collapse because of shear effects. This mass function for 
(Tc = 0.74 and 6c = 1.686 is plotted in figure (6) as a long- 
dashed line. The number of objects with A < 1 has again 
greatly decreased compared to the second axis collapse. 
As could be expected, since now the shear always slows 
down the collapse, the mass function is now much below 
the standard PS prediction. 

In the next section we investigate the generic behavior 
of the mass function obtained in this section (Eq. [27]) in 
the low-mass regime. 

3.2. Low Mass Behavior of the Mass Functions 

In the previous section, we postulated a specific form for 
the collapse epoch of a given fluid element, as a function 
of 6 and the shear along the third axis. In this section, we 
study in a more general way the influence of the shear on 
the formation of low mass objects. We consider that the 
function ac has the following asymptotic behavior when 
X ^ ac 

ac — > (o-c - (31) 

where e > 0. This parameterization is very general, and 
makes no assumptions on the general shape of the func- 
tion dc- We calculate analytically the low-mass (high A) 
behavior of the mass function as a function of e. After 
some algebra (cf. Appendix B), we obtain two character- 
istic regimes, valid in the limit A 1 

^ocA-i"^) iie>^ (32) 
$ oc A"^ if e ^ ^ (33) 
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Note that for the case e = 1, the low-mass behavior is the 

same as for the standard Press & Schechter mass function. 
The domain e > 1 corresponds to shaUower slopes than the 
standard PS mass function. However, the Zel'dovich ap- 
proximation corresponds to e = 1. Therefore, if one con- 
siders that this first order approximation overestimates 
the collapse time, this imposes that 

< e < 1 (34) 

We plot in figure (7) the mass functions we obtained 
using the collapse epoch from formula (24) for different 
values of e < 1. We set the parameter = 3 for each 
curve in figure (7). This parameter determines the point 
where the change of slope occurs in the mass function (in 
our case, this corresponds to A ~ 5). 

This result is extremely important as it shows that 
shear effects can easily modify the low-mass behavior of 
the mass function. The regime e < 1/6 is of primary im- 
portance, because it is independent of the specific form 
of the function ac, which is unfortunately an unknown 
quantity. In this case, for a power law power spectrum, 
one obtains in the low-mass end of the mass function 
$(M) DC M"'+^ where n is the index of the power spec- 
trum. Assuming a constant mass-to-light ratio, this leads 
to a faint end slope of the luminosity function a = n + 1. 
Note that the corresponding power law for the standard 
PS mass function would have been a = — 2 + (n + 3)/6. 
The observed luminosity function of field galaxies exhibits 
this power law behavior at low-luminosity with a ~ — 1 
(Loveday et al. 1994; Efstathiou, Ellis & Peterson 1988). 
It has been previously argued (Blanchard et al. 1992) that 
for a CDM power spectrum with n ~ — 2 at these scales, 
it is impossible to reconcile the PS formalism (a ~ —1.83) 
with the observed luminosity function without imposing 
some ad-hoc procedures, such as a variable mass-to-light 
ratio, or some other biasing mechanism due to the dissipa- 
tive physics of the baryonic gas. Here, we recover a ~ — 1 
with a power spectrum n ~ — 2, using only gravitational 
processes, and more precisely by considering the effect of 
the shear. We claim here that the a CDM-like power spec- 
trum is able to reproduce the low-mass end of the lumi- 
nosity function, when the shear and tidal effects are cor- 
rectly taken into account. Note that the PS formalism is 
better suited for field galaxies than for galaxies in clusters 
where tidal effects are much stronger. With this formal- 
ism, the clusters are considered as single objects whose 
internal structure can not be described. 

4. Conclusion 

In this paper, we extend the PS formalism to more com- 
plex, and perhaps more realistic dynamical models than 
the standard spherical model. We study the infiuence of 
the shear and tidal effects on the mass function of cosmic 
structures. 
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Fig. 7. Mass function obtained for the third axis collapse 
and for e=l, 1/2, 1/4, 1/6 and 1/100 (solid, dotted, dashed, 
long-dashed and dotted-dashed lines respectively). Note that 
for e < 1/6 the slope at low mass remains equal to —7. 

The Lagrangian evolution of a given halo from the ini- 
tial volume up to the final dense and virialized object, 
can be divided into three stages: the collapse of the first, 
second and third axes. Each stage corresponds to a cer- 
tain geometry of the halo formed (sheet-like, filamentary 
and quasi-spherical) and also to a certain collapse epoch. 
In the Lagrangian formalism, this collapse epoch depends 
purely on the initial shear and density contrast. The mass 
function of these halos is then obtained from the collapse 
epoch and from the statistics of the initial parameters. 

We investigate in this paper the influence of the initial 
shear on the dynamics of the forming halos. In the case of 
sheet-like objects, the Lagrangian theory of gravitational 
dynamics allows us to compute the collapse epoch of the 
first principal axis as a function of 5 and 0-3 and 0-2, the 
initial density contrast and two eigenvalues of the initial 
shear tensor. We use for that purpose a new approximation 
presented by Audit & Alimi (1996). In this case we obtain 
a mass function which is well fitted by a PS mass function, 
with 5c = 1.4. This result is explained by the fact that 
shear and tidal effects always accelerate the collapse of 
a fluid element along its first principal axes. We therefore 
obtain more high-mass objects than the standard PS mass 
function (see also Monaco 1995). 

We then investigate the mass function of objects that 
collapse along their second and third principal axis. We 
considere simple and general models since Lagrangian the- 
ory is unable to describe the evolution of a fluid element 
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beyond the collapse of the first axis. We find that the 

resulting mass function in the case of the second axis col- 
lapse is very similar to the standard PS mass function 
{6c ^ 1.686). 

Finally, wc derive analytically a new formula for the 
mass function of objects resulting from the collapse of 
their third axis. This depends on two parameters, namely 
a density threshold Sc and a shear threshold CTc- We recover 
the standard PS formula in the case ac +oo. For val- 
ues of fJc comparable to unity, the peak of the multiplicity 
function is lower than the peak obtained by the standard 
PS approach. This could explain why such an effect has 
been detected in numerical simulations (Efstathiou et al. 
1988). We also derive analytically the low mass regime of 
the mass function computed in a wide class of models. We 
show that it is possible to obtain a low-mass slope very 
close to the faint-end slope of the luminosity function of 
field galaxies. 

In all the dynamical models we use, the resulting mass 

functions show common features such as a high-mass ex- 
ponential cut-off and a low-mass power law. The exact 
shape of the multiplicity function around the peak is, how- 
ever, directly related to the underlying dynamical model. 
This has to be tested in numerical simulations, where the 
mass function should correspond to a dynamical model 
that drives the collapse of every halo (at least in a statis- 
tical sense). We postpone this numerical work to a com- 
panion paper. 
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A. Collapse Epoch using our Dynamical Model 

In this appendix we present the collapse time obtained 
from the AA dynamics. The following formulae were ob- 
tained by fitting a thousand points for 5 > and ^ < 0. 
The fits are accurate to about 1%. 
For 6 > 0, we obtain 



x'^ - 1 



where x G [1,2] is defined by 




The functions /+ and /~ depend only on a^/S, and are 
proportional to 5, according to the scaling property de- 
rived in section 2. 



f+ = d 



0.2(1-5+)) 



g+ = {l- Was/d)-'/^ 
f- = 6 0.81^-0.18(1 -g-) 

g- = {l-4.5a3/Sr' 
Note that the spherical model is recovered for vanishing 
shear. The Zel'dovich solution, which is exact when (T3 = 
—25/3 and c72 = —(73/2, is also recovered. 
For 5 <0, we have 



..So 



1 



8-a;' 



,3o 



T -l + ^c-g„_^/'''+g„_^/ 
where x G [1,2] is the same as for 5 > and a is given by 



1 



0.5 



1/2 



3(73/5-1, 

The functions /+ and are now given by 

5+ = (1 + 15 ((73/5 -1/3))-^/^ 

- i 
6 3, 

.g- = (1 + 4.5 ((73/(5 -1/3))-^ 
When the shear is much larger than the density contrast 
both T+ and T~ converge toward the same value as ex- 
pected. 



r 



6 0.81 



0.18(1-5") 
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B. Computation of the low-mass behavior of the 
mass function for the third axis collapse 

We present in this appendix the derivation of the low-mass 
behavior presented in section (3.2). From equations (20) 
and (26) the universal mass function can be written as 

p+oo 

$(A) oc / e-" /^P^,iudciuA)ydc{uA)diy 

JSJA (Bl) 

We want to study the behavior of the previous integral 
for the case where A — > oo. The collapse time has the 
behavior given by equation (31) which allows to calculate 
the following limits 



d'c{x) — > CTc and (Jc{x) — > -x (^+^)/^ 



OO £ 



(B2) 

We then divide into two the integral giving the mass 
function 

Too = IT 1 -'"^^-P,^ (^A))^.e^j (^A) \u (B3) 
T^=Q^]-'^/''V,^i'^a~i{iyA))^^dii,.A)\,. (B4) 

Using the limits given above, we have for large values of 
A 

p-\-oo 

Ioo~ / l-^'/^P^^(i/aj)i.e2^(j.A)-(^+-)/^ri/ 

^ (B5) 

Therefore, the asymptotic behavior of I^o when A — > oo 
is 

~ A-^ (B6) 

Changing back from ly to S and using the limits given in 
equations (B2), can be written as 

' oc \ 

Y.fn]d5 (B8) 



L. 



A 

So 



\n=2 



where /„ = a„ ^ ^^2^+3 and a„ are numerical constants. 
[P^^ and the exponential have been developed in series). 
If e > l/2(n + 1), / fndS diverges and we have 

fnd5 ^A-'-^ (B9) 

lie < l/2(n-|-l) < 1/6, / /„(i(5 converges and therefore 

i-A 

fndS ~ A-2"+3 < A-^ (BIO) 

Keeping only the leading terms from equations (B6), 
(B9) and (BIO) gives the asymptotic behavior of ^>(A) 
(Eqs. (32) and (33)). 
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